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ABSTRACT 

Aims. We quantify the intrinsic width of the red giant branches of three massive globular clusters in M3 1 in a search for metallicity 
spreads within these objects. 

Methods. We present HST/ACS observations of three massive clusters in M31, G78, G213, and G280. A thorough description of the 
photometry extraction and calibration is presented. After derivation of the color-magnitude diagrams, we quantify the intrinsic width 
of the red giant branch of each cluster. 

Results. This width translates into a metallicity dispersion that indicates a complex star formation history for this type of system. For 
G78, <r [Fe/H ] = 0.86 ± 0.37; for G213, 0.89 ± 0.20; and for G280, 1.03 ± 0.26. We find that the metallicity dispersion of the clusters 
does not scale with mean metallicity. We also find no trend with the cluster mass. We discuss some possible formation scenarios that 
would explain our results. 

Key words, galaxies: evolution - Local Group - galaxies: star clusters - globular clusters: general - globular clusters: individual 
(G78, G213, G280) 



1. Introduction the halo of M31 constitutes the only evidence for a metallicity 

spread, as derived f rom the width of its red giant branch (RGB) 
_CS ; There are now a few dozen Galactic globular clusters for which (Mev lanet al200ih . The parallel with co Cen is immediate. Both 
high resolution spectroscopy has allowed insight into their dusters exWbit a metal abundance rang e of about 1 dex, Gl be- 
star-to-star chemical variations ( jGratton et al. | |2004j). Omega ing twice as massive as u Cen The two clusters differ in a num . 
Centaun (co Cen) appears to be the only case where such vana- ber of ways thoughi which are nkely important clues to the un- 
tions are clear for all elements. Several scenarios have been pro- derstanding of their formation. For example, Gl and co Cen have 
posed to explain these observations: The high primordial mass very differe nt locations with respect to their host galaxy - Gl is 
of co > Cen might be enough to induce secondary star forma- a halo cluster located in pro j ec tion at ~ 40 kpc from the center 
tion ( Dopita&Srnith| | 1986|) . Models of chemical evolution have of M3 j (nearly the distance bet ween our Galaxy and the LMC), 
considered self ennchment_scenarios d lkuta & Anmoto| |200fJ while w C en is now at ~ 6 kpc from the Galactic center and at ~ 
iTsuiimoto & Shigevamal 120031) . Primord ial chemical inhomo- j f rom me G a i ac tic plane 
geneities have also been investigated (e.g. lKrafll (IT994)). Finally, 

the possibility that co Cen is the nucleus of a now-dissolved nu- Although co Centauri is by far the brightest and most mas- 

cleated dwarf galaxy has been raised by a number of studie s srve globular cluster in our Galaxy, G 1 may not be the only such 

dZinnecker e t all l 19881 : iLee et alj [l999l: iHilker & Richtleiil2000l) . massive cluster belonging to M3 1 . There are at least three other 

Our knowledge of extragalactic globular clusters is currently bright clusters w hich have central velo city dispersions larger 

far less advanced. At the moment, the globular cluster Gl in than 20 km s -1 dDiorgovski et all Il997b . They offer a unique 

opportunity to investigate the degree to which the mass of a 

Send offprint requests to: I. Fuentes-Carrera globular cluster influences its chemical evolution. In the veloc- 

* Individual photometric measurements will be made available at the ity dispersion versus absolute magnitude relation (the equivalent 

CDS. for globular clusters of the Faber-Jackson relation), the two se- 
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quences of Galactic and M31 systems show no detectable dif- 
ference in the slope and zero-point of their correlations. Thanks 
to the similarity between the Galactic and M3 1 globular cluster 
sequences, any conclusion drawn from additional M3 1 clusters 
may represent general (universal) characteristics. 

We present here the analysis of images obtained with 
the Advanced Camera for Surveys aboard the Hubble Space 
Telescope of those three massive and bright clusters in M3 1 . Our 
goal is to quantify the intrinsic width of their red giant branches 
in order to search for any metallicity spread. 



2. Observations and Data Reduction 

2.1. Sample selection 

Velocity dispersions are a good indicator of dynamical mass 
for old stellar systems with nearly constant mass-to-light (M/L) 
ratios, such as globu lar clusters dDiorgovski & Mevlanl 1 1994; 
Diorgovs ki et alll997l) . Therefore, our selection criterion for this 
work operates via the cluster total luminosities and velocity dis- 
persions. 

As mentioned previously, apart fro m Gl, three other clusters 
among the 21 objects in the sample of Diorgovski et~al i 19971) 
have central velocity dispersions, cr, greater than 20 km s~ l . 
Figure Q] shows the location of the three clusters in M3 1 . All 
three of them are located off the major axis of M31 towards 
the minor axis. G78 is located (in projection) on the edge of 
a dust cloud (object D48 in Hodge 1981) which lies northwest 
of the cluster. G213 coincides in projection with one of M31's 
spiral arms but does not lie near any particular dust cloud or 
stellar association. G280 is also located at the edge of a dust 
cloud (object D642) lying to the northeast of the cluster. Table 
Q]lists some of the properties of these clusters. Columns (1) and 
(2) give the coordin ates of each cluster, column (3 ) the magni- 
tude in the V filter dHuchra. Brodie & Kentlll99ll) . column (4) 
the projecte d distance from the center of M31, column (5) the 
metallicity ( Huchr a. Brodie & Kent 199 lb . column (6) the ve- 
locity dispersion dDiorgovski et all 1997b and column (7) the red- 
dening value (Jablonka, Alloin & Bica 1992 for G78; this work 
for G213; Frogel, Persson & Cohen for G280). 



2.2. Observations 

Observations of the three clusters and their surrounding fields 
were obtained with the High Resolution Channel (HRC) of the 
Advanced Camera for Surveys (ACS) aboard HST (Cycle 12; 
program ID 9719). The HRC has a plate scale of 0.027" pixel- 1 
and a field of view (FoV) of 26" x 29". Four images of each 
globular cluster were taken using the F606W (close to V) and 
F814W (close to I) filters with total integration time of 505 x 4 
= 2020 s and 715x4 = 2860 s, respectively. In order to be able 
to detect cosmic rays and other artifacts as well as to sample the 
PSF optimally, images were dithered following the four-point 
dither pattern known as ACS-HRC-DITHER-BOX with a step 
of 0.15" eachfl 

For each cluster, we used the FLT images downloaded from 
the HST archive. The instrumental signature has been removed 
from these images (bias, dark and flat-field correction) by the 
on-the-fly-reprocessing (OTFR) performed by CALACS. These 
FLT images were then corrected for geometric distortion and 



1 See the ACS Phase II Proposal Instructions 
http://www.stsci.edu/hst/acs/proposing/dither for further details 



combined using the Mul tiDrizzle task of the S TSci pipeline in- 
corporated in PyRAF ( Koekemoer et al 2002). This procedure 
eliminated all of the cosmic rays and bad pixels on the images. 
We have used the default values for all of the parameters in 
the MultiDrizzle task. We used "median" as the IMCOMBINE 
method in order to add up the four images per filter. Due to a gra- 
dient in the background of the ACS field, images were not sky- 
subtracted. Background variation was dealt with during the pho- 
tometry extraction using DAOPHOT (Section [3]). No oversam- 
pling was done when drizzling the images into the final image to 
avoid introducing a second redistribution of noise in the resam- 
pled pixels. Figure|2]shows the final MultiDrizzled image of each 
cluster in the F8 14W filter. The long, dark feature on each image 
is the shadow of the occulting finger of the ACS/HRC corona- 
graph. A quick comparison of these images shows that G78 is 
bigger than G213 and G280 which are rather similar in size and 
structure. All of them have a dense core, though proportionally, 
the core of G78 is smaller with respect to the cluster as a whole. 
The field surrounding G213 is the most populated. Both of the 
fields surrounding G78 and G280 are sparse, although that of 
G78 seems to be somewhat less populated. 

3. Photometry 

3.1. Photometry Extraction 

Photomet ry was determ ined using the DAOPHOT/ ALLSTAR 
routines (Stetson 1994). Since MultiDrizzled images are pro- 
duced in units of electrons second^ 1 pixeV x , they need to be 
multiplied by the exposure time and divided by the gain factor to 
get back to the original ADU pixeV 1 units used by DAOPHOT 
and ALLSTAR. In order to do so we need to consider how many 
original images are included in the drizzled images. If N images 
have been combined by the drizzle routine, then the read noise is 
reduced by a factor between 1 and "V7Y, and the gain is increased 
by a factor between 1 and N. Unfortunately these numbers are 
not necessarily constant because any one pixel in the drizzled 
image can come from as few as 1 or as many as N original 
images. One must also consider the fact that the MultiDrizzle 
task redistributes the noise in each pixel of the image, while 
DAOPHOT assumes that the noise properties of an image are ap- 
proximately constant everywhere. In order to take these effects 
into account in the derivation of precision photometry, we fol- 
lowed P. Stetson's suggestion and modified our initial values for 
the READ OUT NOISE and the GAIN in both the DAOPHOT 
and ALLSTAR routines in order to have x 1 values close to 1 . 
Object detection was done with DAOFIND on the combined 
MultiDrizzled image. Point-spread functions (PSFs) were con- 
structed with DAOPHOT and stars were measured using the 
ALLSTAR routine. The PSF of each cluster on each filter was 
built from bright isolated stars in the FoV. For G213 and G280, 
the PSF for the F606W image was best determined using a few, 
bright and isolated stars (23 for G213 and 25 stars for G280) 
while for the F814W image, the PSF was better computed using 
more stars (52 for G213 and 47 stars for G280). For G78, both 
PSFs were best determined using a small number of stars (22 
stars for both images). For the three clusters, the stars used to es- 
timate the PSF cover the whole FoV except for the central parts 
of the FoV (about 2 arcsec from the center of the FoV). As a first 
approach, we used a PSF that varied linearly across the FoV in 
case there still remained a slight variation of the pixel scales due 
to the optics of the ACS. Through the analysis of residuals after 
star subtraction, we found that the photometry was best achieved 
using a spatially constant PSF, which supports the fact that the 
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Fig. 1. Location of globular clusters in M31. Dark purple triangle corresponds to the location of G78, light blue circle corresponds 
to G213 and dark r ed square sho ws the location of G280. Numbered squared regions correspond to charts in the Atlas of the 
Andromeda Galaxy (iHodgdl 19811) 

Table 1. Properties of observed clusters. 



Cluster 


a (2000) 


8 (2000) 


V(mag) 


DM3lC) 


[FelH\ 


<r(km/s) 


E(B - V) 




(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


G78 


00 41 01.2 


+41 13 45.5 


14.2 


19.2 


-0.92 


25.48 


0.23 


G213 


00 43 14.4 


+41 07 20.5 


14.5 


11.2 


-1.08 


20.50 


0.10 


G280 


00 44 29.5 


+41 21 35.8 


14.3 


20.5 


-0.70 


25.94 


0.10 



References: (3) and (5) Huchra, Brodie & Kent 1991; (6) Djorgovski et al. 1997; (7) Jablonka, Alloin & Bica 1992 for G78, this work for 
G213, Frogel, Persson & Cohen for G280. 



Fig. 2. MultiDrizzled images of G78 {left), G213 (middle) and G280 {right) and surrounding field stars in the F814W filter. North is 
to the top, East is to the left. The dark feature seen on the northeast corresponds to the occulting finger of the ACS. Field of view 
(FoV) for each cluster is 17.3"x 17.3 ". 
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pixel variation was properly corrected by the MultiDrizzle pro- 
cedure. 



3.2. Photometric Calibration 

Once the instrumental photometry was extracted, it was n eces- 
sary to convert it to the STMAG system dKoornef et allll986l) . In 
order to do so, the MultiDrizzled nature of the images must be 
taken into account. Our MultiDrizzled images were obtained by 
drizzling N FLT images into a single image. FLT images have 
units of electrons while the MultiDrizzled image is in units of 
electrons s . Photometry was extracted using the DAOPHOT 
and ALLSTAR packages which work with ADU . As previously 
mentioned, in order to work with a MultiDrizzled image using 
these packages, we multiplied the original MultiDrizzled image 
by the exposure time of one FLT image, t exp , and divided it by 
the gain of the CCD, gain cc ^ We shall call the resulting image, 
the modified Multidrizzled image. 

3.2.1. Charge Efficiency Transfer corrections 

Photometric losses due to the problem of charge transfer effi- 
ciency (CTE) of the CCD need to be quantified before convert- 
ing instrumental magnitudes to any photometric standard sys- 
tem. For the ACS/HRC, the CTE correction is given by the for- 
mula in the ACS Data Handboolfl In order to have the cor- 
rect input values, the SKY brightness obtained with ALLSTAR 
from the modified MultiDrizzled image in ADU, must be mul- 
tiplied by the gain of the CCD. On the other hand, the FLUX 
of the star (in electrons) is the original flux on the FLT image. 
This flux is equal to the flux associated to the magnitude given 
by DAOPHOT, MAG modDRZ (star), times the gain of the CCD, 
gain ca i. Considering that DAOPHOT adds a zeropoint of 25.0 to 
all magnitudes, the resulting flux is given by 



FLUX (electrons) = io-(MAC m „ JDRZ (s,ar)-25.0)/2.5 x (1) 

Once the SKY and FLUX values have been corrected, the 
resulting CTE correction (YCTE) will be that for a single FLT 
image. 

3.2.2. Aperture Corrections 

Our aim is to transform the magnitudes obtained from 
ALLSTAR, which we will call PSF magnitudes, into magni- 
tudes in the STMAG system. In order to do so, several aperture 
corrections need to be made. Before being converted into mag- 
nitudes in an aperture with "infinite" radius, PSF magnitudes 
need to be corrected to magnitudes in an aperture of radius 0.5" 
(18.5 HRC pixels). Since the aperture corrections used to trans- 
form instrumental magnitudes into the STMAG system have 
been computed for single FLT images (in units of electrons), we 
need to "divide" all PSF magnitudes derived from the modified 
MultiDrizzled image, MAG moc {DRz{PS F) in ADU, by the gain of 
the CCD. 



MAG fit (PS F) — MAG mo dDRz(PS F) - 2.5 log(gain ccd ) (2) 



We then select 10-20 uncrowded and bright (F606W and 
F814W magnitudes brighter than 16.5 mag and 15.5 mag, re- 
spectively) stars on the MultiDrizzled image and measure their 
total magnitudes inside a radius of 18.5 pix, MAG mo dDRz(^.5), 
using PHOT in DAOPHOT. Since this magnitude value will be 
that corresponding to the modified MultiDrizzle image (whose 
flux is in ADU), we need to convert it into magnitudes on the 
FLT image (whose flux is in electrons). 



MAG FLT (18.5) = MAG 



(18.5) - 2.5 \og(gain ccd ) (3) 



modDRZ 1 

For each image, we calculate the mean value of the differ- 
ence between these two magnitudes. The mean value of this dif- 
ference is the aperture correction, AP_CORR(18.5), for each FLT 
image. This value is then added to all of the PSF magnitudes 
of the FLT image, MAGflt(PS F). The resulting magnitudes, 
MAG apcorr , need to be converted from this "intermediate" aper- 
ture to "infinity" or observ ed magnitude, OBMAGFLTmf- This is 
done using equation (2) in ISirianni et all (120051) . Strictly speak- 
ing, the Sirianni et al. corrections to infinity are only applicable 
to the drizzle d images, not the FLT images. Nevertheless, in a 
recent paper, ISaraiedini et all (120061) show that there is no sig- 
nificant difference between them. In order to use this equation, 
MAG a p Corr needs to be set to COUNT JiATE in electrons s~ l 
through the following equation 



COUNT _RATE(electrons S -1 ) 



10 M/lG„„ ( . o , T (18.5)/2.5 



(4) 



Once all magnitudes have been set to the observed mag- 
nitude, OBMAGFLTinf, they need to be corrected by the CTE 
correction previously computed. The OBMAG f ma i will be used 
to transform the instrumental ma gnitude system to th e STMAG 
system following equation (7) of lSirianni et all (120051) . 



3.3. Conversion to STMAG absolute magnitudes 

Finally the observed STMAG magnitudes need to be converted 
into absolute STMAG magnitudes taking into account the dis- 
tance at which the clusters lie and thei r reddening v alue . This 
was done following the proc edure by Brown et all ( 2005b . A 
value of (m - M) = 24.43 dFreedma n & Madore 1990) was 
used for the true distance modulus. For each cluster, two redden- 
ing values were considered: 0. 1 and 0.23 The value E(B-V)-0. 1 
is the Milky Way line of sight reddenin g in the direction of M3 1 . 
It is th e value estimated for G280 by iFrogel. Persson & Cohenl 
(119801). E(B - V)=0.23 is th e value estimated for G78 by 
iJablonka. Alloin & Bical d!992l) . There is no reddening value in 
the literature for G213. 



2 This can only be done if ALL FLT images used in the drizzling 
process have the same exposure time, t exp 



3 chapter 6, section 1.5 



4. Statistical Selection of Stars 

Since we are interested in determining the width of the RGB 
of each cluster, we need to eliminate any spurious photometric 
detections that might contaminate our results. In order to do so, 
stars were statistically selected to eliminate objects with dubious 
values for the error, x 2 and sharpness. This was done prior to the 
photometric calibration, that is with the star magnitudes given in 
instrumental magnitudes.Figures[3]|4]and|5]show the output of 
DAOPHOT for the total ACS field of the three clusters. For all 
three parameters, most points seem to follow a global trend. For 
instance, in the case of G280, error values increase from ~ 0.02 
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at 16.0 F814W to 0.4 at 22.0 F814W. There is however, a per- 
centage of outliers for each distribution for all magnitude values. 
We considered these points as possible spurious photometric de- 
tections to be removed before any analysis. A statistical selec- 
tion was done on the estimate of the dispersion and the mean 
value of the distribution of each parameter. The mean value of 
the distribution was computed recursively; after each computa- 
tion, stars beyond 3 sigma of the distribution were eliminated 
before the next computation of the mean. A bin of 0.7 mag was 
used for parameters in the F606W filter and a bin of 0.5 mag in 
the F814W filter. The computation of the mean was performed 
until the difference between subsequent values of the mean value 
of each parameter became lower than 1% of the value. Once 
this final mean value was determined, we considered the highest 
value of the dispersion of each parameter (since the computation 
was done independently for each magnitude bin the dispersion 
values vary from one bin to the other). This value for the disper- 
sion was multiplied by a factor between 1.0 and 3.5 in order to 
delimit the edges of the distribution. These factors were chosen 
in order to draw the boundaries that follow the bulk of the dis- 
tribution of each parameter value. The solid lines in Figures [3] [4] 
and [5] show the upper limit for the statistical selection of stars; 
the lower limit for the selection was also considered in the case 
of the sharpness. The dashed lines show the average value of 
each parameter. For all clusters, stars whose x 1 values were over 
3.5 for both filters were discarded before any statistical selec- 
tion was done, in order to facilitate the statistical selection and 
to eliminate stars that had undetermined sky flux to be used for 
the CTE correction. 

Following this selection, we were left with 7770 stars from 
the initial 8139 stars for G78. For G213, 13275 stars from 14781 
initial stars were selected. For G280, 11204 stars from 12409 
were kept. 

5. Artificial Stars Experiments 

To evaluate the photometry extraction and to quantify the effects 
of stellar crowding, we performed ten sets of artificial star exper- 
iments per cluster. These were done using instrumental magni- 
tudes. In order to do so, we constructed the CMD of each cluster 
in instrumental magnitudes. Left panels of Figure [6] show the 
CMD derived for each cluster, right panels show the CMD built 
considering all the stars in the ACS FoV -after statistical selec- 
tion. For the cluster CMDs, stars were selected within a circle 
centered on each cluster. This was done through the analysis of 
the radial distribution of stars on each FoV. Figure [7] shows the 
radial density profile for the entire FoV centered on each cluster. 
The outer radius of the cluster, R oul , was chosen to be the radius 
at which the density profile dropped significantly before reach- 
ing a fairly constant value in order to minimize contamination 
from field stars. This plateau is assumed to represent the density 
of the surrounding field. Values for R out are shown in Table [4] 
We then derived the fiducial RGB for each cluster by fitting the 
mean locii of stars between 15.8 and 18.2 mag in the F814W 
filter. The faintest limit of the magnitude interval was chosen to 
avoid the mixed stellar population of the red clump. 

For each experiment, we selected 384 stars along the fiducial 
RGB of each cluster from 15.8 to 18.2 in F814W in magnitude 
bins of 0.05 to fully cover the RGB interval. This resulted in 8 
stars per magnitude bin per experiment and 80 stars per magni- 
tude bin for the ten experiments as a whole - in total, 9216 artifi- 
cial stars were used for the experiments. The artificial stars were 
added to both the F606W and F814W HRC/ACS MultiDrizzled 
frames with the constraint that no two artificial stars be within 
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Fig. 6. Left panels: Color-magnitude diagrams (CMD) for G78 
(top), G213 (middle) and G280 (bottom) considering all stars 
within R ou , -values given in Table|2] Right panels: CMDs consid- 
ering all the stars in the field of view (FoV) of the MultiDrizzled 
HST-ACS images of G78 (top), G213 (middle) and G280 (bot- 
tom). 



1.5 PSF radii of each other and avoiding the edges of the im- 
ages. The number of artificial stars considered per experiment 
was chosen to be small enough to avoid introducing important 
density variations in the FoV of the clusters. The number of ar- 
tificial stars introduced in the FoV is at most 5% of the total 
number of real stars in the field. The photometry from the resul- 
tant images was computed following the exact same steps used 
to derive the photometry of the observed images (see section 
13.11 ). The percentage of recovered stars for all ten experiments 
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14 16 18 20 
F606W 




Fig. 3. DAOPHOT/ALLSTAR output for G78. Top: Error vs. magnitude for both F606W and F814W instrumental magnitudes. 
Middle:x 2 vs. magnitude. Bottom: Sharpness vs. magnitude. Solid lines show the upper limits for stars used to build the CMDs of 
the cluster and the field; dashed line shows the averaged value after iteration. In the case of the sharpness a lower limit was also 
used to select stars. 
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Fig. 4. DAOPHOT/ALLSTAR output for G213. Same as Figure[3] 



14 16 18 20 
F606W 
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per cluster is 97.94 %, 98.13 % and 98.29 % for G78, G213 
and G280, respectively. The missing fraction are artificial stars 
whose photometry could not be extracted with DAOPHOT be- 
cause they are located too close (< 3px) to real stars in the FoV. 

Figure [8] displays the input artificial stars along the fiducial 
RGB and the recovered RGB stars for all the experiments in each 



cluster. No statistical selection of output stars was done for this 
comparison. For the three clusters, the input sequence is well re- 
covered. Less than 5.0 % of recovered stars fall below the magni- 
tude 18.2 in F814W -for G78 it is less than 2.5 %. For this same 
cluster, less than 1.0 % of the recovered stars fall to the right 
(and/or left) of the fiducial RGB - between 1 and 2 magnitudes 
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Fig. 5. DAOPHOT/ALLSTAR output for G280. Same as Figureg] 



in color. The same occurs for the recovered stars of G280. For 
G213, there is a higher percentage of stars (~ 3.0%) falling out- 
side the "main" body of the RGB than for G78 and G280. This is 
probably due to the fact that the FoV of this cluster is more pop- 
ulated than those of G78 and G280. In all cases, the number of 
outliers from the recovered RGB is small. These are stars whose 
recovered magnitudes differ significantly from the input mag- 
nitudes because they happen to overlap with the original (real) 
stars in the FoV. Although the artificial stars were placed so that 
no two artificial stars are within 1.5 PSF radii of each other, no 
such restriction was imposed for stars already in the FoV. 

An evaluation of the input minus the recovered magnitudes 
in the F606W and F814W filters according to the position of the 
star from the center of the cluster showed that this difference in- 
creases at small radii (R ~ 4.8", 3.5" and 4.0" for G78, G213 
and G280, respectively) introducing an artificial brightening for 
stars near the center of each cluster. Figure [9] shows the differ- 
ence between the input magnitudes of the artificial stars and their 
recovered magnitudes (magJLiff) according to their radial posi- 
tion with respect to the center of each cluster. An important de- 
parture from zero is seen for low values of R. For higher values, 
some stars do "migrate" in magnitude, but they are few if com- 
pared to the bulk of the stars which fall around zero. In order to 
evaluate this brightening for the stars in each cluster, we com- 
puted the average value of magjliff and its dispersion within 
R out . We considered all stars within an annulus of outer radius 
R oul and whose inner radius we varied until the absolute aver- 
age value of mag-diff was less than 0.02 mag and the disper- 
sion value was close to 0.05 mag or less. Results are shown in 
Table[2] This final radius was set to be the inner radius R,„ of the 
cluster. All stars lying within this radius will not be taken into 
account in the analysis that follows. 

Sources of error, e.g., such as residual flat-field non- 
uniformities and residual dark current, are not included in our 
photometric error estimates. However, those should not dramati- 
cally increase our estimates. There is also a possibility that small 



Table 2. Average value and dispersion of the difference between 
the input magnitude and recovered magnitude between specific 
radial intervals for each cluster 



Cluster 


radial interval (") 


mean magnitude difference 


dispersion 


G78 


5.1-6.4 


-0.0203 


0.0454 


G213 


3.7-5.3 


0.0343 


0.0551 


G280 


4.2 - 5.6 


0.0425 


0.0528 



differences exist between the PSF model used in the artificial star 
experiments and the profiles of the real stars. This source of er- 
ror is difficult to quantify; however, we believe it is likely to be 
negligible because of the spatial resolution of the ACS/HRC in- 
strument which provides good sampling of the PSF (FWHM ~ 
3 pixels). 



6. Luminosity function 

The luminosity function (LF) for each cluster and its surround- 
ing field stars was derived using a bin of 0.15 in magnitude^. 
Field stars were selected within a concentric annulus with an in- 
ner radius 1 .0" larger than the outer radius of the cluster, R out 
-defined in section [5] This was done to reduce the number of 
cluster stars that could be contaminating the field. The outer ra- 
dius of the ring containing field stars was chosen 3.0" larger than 
its inner radius. The LF of each cluster was computed consider- 
ing all of the stars between the inner radius, Ri„ and outer radius, 
R out defined in the previous section. Figure [10] shows the lumi- 
nosity function for each cluster and its corresponding surround- 
ing field. Figure [TT] shows the resulting luminosity function of 
each cluster subtracted from the contribution of the stars in the 
surrounding field. For G 78, an E(B - V) of 0.23 was c onsid- 
ered initially (taken from Jablonka, All oin & Bical (119921) ). This 



4 Absolute STMAG magnitudes are considered from this part of the 
analysis onwards. 
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G78 



Table 3. Magnitude and color values for the boundaries of the 
intervals considered to fit the RGB of each cluster 



10 



radius (arcsec) 



G213 



radius (arcsec) 



5 10 

radius (arcsec) 

Fig. 7. Density profile for all stars in the FoV of G78 (top), G213 
(middle) and G280 (bottom) after statistical selection (see sec- 
tion [4}. Dashed lines delimit the region used to construct the 
CMD of each cluster. Dot-dashed lines delimit the region of the 
FoV used for the CMD of the surrounding stars (see section l7T2l i. 
Apparent density holes near the center of the clusters are due to 
star crowding. 



Cluster 


lower mag 


upper mag 


lower color 


upper color 


G78 


-2.2 


0.685 


none 


0.8 


G213 


-2.2 


0.685 


-0.7 


0.5 


G280 


-1.8 


0.685 


none 


0.6 



7. RGB width determination 

7.1. Final CMD 

The final CMD of each cluster was built considering only statis- 
tically selected stars within the ring between R,„ and R out . Figure 
[T2l shows the superp osition of the stand ard RGB sequences and 
HB loci derived by iBrown et all d2005l) on the CMD of G78, 
G213 and G280. The straight dotted line shows the position of 
the red clump of each cluster. For the three cluste rs, the red 
clump falls at the same position as the HB loci of IBrown et all 
(120051) . 

At this point, the CMD of each cluster should be largely free 
of spurious detections. Nevertheless it will still be affected by 
the presence of field stars falling within the area of each cluster. 
In order to estimate the contribution of field stars to the RGB 
of each cluster, we also constructed the CMD of the stars in the 
surrounding fields. These stars were chosen within the outer an- 
nulii described in section[6] The left panels of figure[T3]show the 
CMDs of G78, G213 and G280. Right panels of the same figures 
show the CMD of their corresponding surrounding fields. 

For the final CMD of each cluster, the fiducial RGB was de- 
rived by fitting the mean locii for stars within certain magnitude 
and color intervals. The upper magnitude limit was set at 0.5 mag 
brighter than the position of the red clump. The lower magnitude 
limit, as well as the color limits, were set in order to get rid of 
obvious outliers of the RGB. Magnitude and color intervals for 
each cluster are shown in Table[3] The fit is shown in Figures fT2l 
and [13] The dark lines in the left panels of these Figures show 
the mean locii of the RGB regions of each cluster. 

7.2. Density distribution of fields stars vs. cluster stars 

In order to determine the actual contribution of field stars to the 
RGB of the cluster, we need to compare the density distribution 
of field stars to that of stars in the cluster. Figure [7] shows the 
density distribution of the clusters and their surrounding field. 
Since the cluster stars have been selected in view of the robust- 
ness of their photometry, the most central stars of the cluster are 
missing. This leaves a hole in the radial distribution of stars. The 
average density of stars in each annulus (associated to each clus- 
ter and to its surrounding field), as well as the surface of the 
cluster are given in Table [4] 



value was then changed to 0.22 in order to match the location 
(in magnitude) of the HB of G78 to the location (in magni- 
tude) of the HBs of G213 and G 280. For G280, E(B - V)=0.1 
(Fro gel. Persson & Cohen|[l980l) . For G213, there are no red- 
dening values in the literature, so the luminosity function was 
computed considering E(B — V)=0.1, which is taken to be the 
standard value along the M31 line of sight. The superposition 
of the CMD of G 213 with the standard RGB of sequences of 
Bro wn et"al d2005) confirms the validity of this choice (middle 
panel of Figure [12}. For each cluster and its surrounding field, 
the red clump lies at ~ 1.1 85. 



7.3. Color difference histograms for clusters and surrounding 
fields 

Figure[T4]shows the color difference histograms around the fidu- 
cial RGB for each cluster and its corresponding field. These his- 
tograms show the difference between the color of a star with a 
certain magnitude and the color value of the fiducial RGB at the 
same magnitude. RGB stars were considered within a magnitude 
interval of 0.5 mag, with its faintest limit set at 0.5 magnitudes 
above the red clump. This value was chosen to avoid the mixed 
stellar population of the red clump, while the brighter limit of 
the magnitude interval (0.5 magnitudes brighter) was chosen to 
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avoid the upper parts of the RGB that could be polluted by AGB 
stars. The magnitude interval for the three clusters goes from 
0.185 to 0.685 mag. Histograms were derived using a color bin 
of 0.085 mag for G78 and G213 and 0.090 mag for G280. The 
size of color bin was chosen to best sample the distribution of 
stars along this part of the RGB according to each cluster. The 
distribution of the field stars was scaled to the same area as the 
cluster stars based on the values of and R ou , shown in Table 
[4] Histograms were obtained by direct bining of stars per unit 
color (straight line histograms in Figure [T4l. So-called "general- 
ized histograms" (filled circles in the same figure) were obtained 
by weighting the contribution of each star to each color bin of the 
histogram with its photometric error, i.e. the value of the photo- 
metric error of each star, determined the weight this star would 
add to the star counts in each color bin of the histogram. For 
all three clusters, the first histogram (simple histogram) and the 
second one (generalized histogram) follow the same trend. Once 
the color distributions of the field stars had been properly scaled 
to the area considered for the determination of the color distri- 
bution of the cluster stars, the field stars were directly subtracted 
from the cluster stars per color and magnitude bin in order to 
derive a color distribution of cluster stars along the RGB. 

8. Analysis of RGB width 

In order to evaluate the contribution of photometric errors to the 
width of each cluster's RGB, we derived the CMD of the recov- 
ered artificial stars after applying a similar statistical selection 
than the one done for the observed stars and described in sec- 
tion |U as well as using the magnitude and color limits shown 
in Table [3] Only artificial stars falling within each cluster annu- 
lus are considered. Figure [15] shows the CMD of these artificial 
stars. 

Considering the same magnitude intervals than the ones used 
in the previous section to compute the color histograms of the 
observed stars, we derived the color histogram of the recovered 
artificial stars. Empty circles in Figure [16] show the generalized 
normalized histogram of the recovered artificial stars per clus- 
ter. The dashed line shows the Gaussian fit. Filled circles in the 
same figures display the generalized normalized histogram for 
the observed stars without the contribution of field stars. Solid 
lines show the Gaussian fits. The resulting histograms for the ob- 
served stars are not centered on zero due to the fact that the fidu- 
cial RGB was derived before the subtraction of field stars, i.e., 
considering both cluster and field stars. Once the subtraction was 
done, the locus of the remaining distribution is slightly shifted 
from zero. For G78, there is a shift of -0.03 mag; for G213, the 
shift equals -0.03 mag; for G280, there is a shift of 0.02 mag. 
Rows 6 and 7 of Table |4] display the values of <x for the two 
Gaussian fits along with their uncertainties. <x ofcs corresponds to 
the distribution of observed stars while cr arl corresponds to that 
of the artificial stars. These values were subtracted in quadra- 
ture in order to analyze the intrinsic width of the observed stars 
distribution for each cluster. The last row of Table [4] show the 
resulting difference, <j int = (cr 2 bs - (r 2 rt ) 1/2 and its uncertainty. 

9. Determination of the color-metallicity relation 
and metallicity dispersion 

Assuming the intrinsic color width of the RGB of these clusters 
is due entirely to a metallicity dispersion, we can estimate the 
value of this dispersion by constructing the color-metall i city re - 
lation using the standard RGB sequences of [Brown et all d2005). 



Table 5. Metallicity dispersion for each cluster 



Cluster ID 


G78 


G213 


G280 


°" [Fe/H] 


0.9419 


0.8955 


1.0337 


error 


+0.3824 


±0.1979 


±0.2587 


0~[Fe/H]/[Fe/H] 


1.0238 


0.8292 


1.4767 


error 


±0.4157 


±0.1832 


±0.3695 



Considering the central magnitude of the interval used to es- 
timate the intrinsic color width of the clusters, we plotted the 
corresponding colo r values for each standard RGB sequence in 
Bro wn et"afl (l2005h versus its metallicity. These points were fit- 
ted using a second order polynomial in order to derive the color- 
metallicity relation for the clusters as shown in Figure [17] The 
resulting color-metallicity relation for the three clusters is: 

[Fe/H] = 0.0952 + 2.8323(V-/) - 12.5153(^-7)^ (1) 

where (V - 7)o is the RGB color at 1=0.435 mag, which is 
0.75 magnitudes brighter than the HB. Using this expression, 
we derived a metallicity dispersion for each cluster. 

First we compute the metallicity for the color (V—I)o, which 
we shall call [Fe/H]o- Then we add to [V - 7)o the intrinsic 
color dispersion computed from the histogram fit in section [8] 
that is cr int , to have (V - 7)o + a- int . This value is inserted into 
equation 1 to derive the "upper limit" of the metallicity disper- 
sion, [Fe/H] up . For the "lower limit" of the metallicity disper- 
sion, we subtract to (V - 7)o the intrinsic color width, that is 
(V - 7)o - cr,„ f , and insert this value into the [Fe/H] vs color 
relation to obtain [Fe/H]i mv . These lower and upper limit are in- 
dicated with crosses in Figure[17] The metallicity dispersion will 
be given by ([Fe/H] up -[Fe/H]i ow )/2 even though [Fe/H]o does 
not lie in the middle of this interval. The metallicity dispersion 
of each cluster is given in Table [5] We give both the resulting 
dispersion and the relative dispersion, o-[f e /H]/[Fe/H] , where 
[Fe/H] is given in Table [TJ 

As mentioned in Section 5, we have performed the photom- 
etry and the artificial stars experiments with great care. The 
latter have been conducted in such a way as to include all 
known sources of photometric error that are readily quantified. 
Nonetheless, there is a small chance that we have underestimated 
the actual photometric errors, in which case, our inferred cluster 
metallicity dispersions would be upper limits of the true values. 

As a comparison, the intrinsic metallicity dispersion of Gl 
(a- = 25.1 km s-1, [Fe/H]=- 1.08) is ±0.5 dex for E(B - V) = 
0.10, and ±0.39 dex for E(B - V) = 0.05 dMevlan et alll2001l) . 
which is equivalent to that of the three clusters within the er- 
ror bars. We note that the G280's metallicity dispersion that 
we derive in this wo rk is much higher than the value found by 
(IStephens et all200ll) in their study involving NICMOS observa- 
tions of metal-rich globular clusters in M3 1 . According to those 
authors, o~[F e /H] ~ 0.2, using the dispersion in the measured 
{J - H) and (J - K) colors. This difference is probably due to the 
fact that for any spread in [Fe/H], the corresponding spread in 
color is twice as small in the infrared (J - K) than in the visible 
(V-I). Also the fact that they work with a smaller sample (~200 
stars versus ~2500 stars per cluster in this work) would statisti- 
cally affect their measurement of any spread in metallicity. 

10. Discussion and Conclusions 

In light of the existing hypotheses made for the case of our 
Galaxy, we now consider some possible explanations for the 
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Table 4. Parameters for the RGB width analysis and sigma values from the Gaussian fits of color distributions of observed and 
artificial stars. 



duster id 








R jn (arcsec) 


5.1 


3.7 


4.2 


R m „(arcsec)v 


6.4 


5.3 


5.6 


Cluster density (stars per arcsec 2 ) 


14 


20 


19 


Field density (stars per arcsec 2 ) 


7 


15 


12 


Cluster surface (arcsec 2 ) 


55 


45 


43 


c obs ± Act obs 


0.1449 ±0.0395 


0.1110 ±0.0261 


0.1439 ±0.0415 


(Tart ± A<r art 


0.0587 ± 0.0142 


0.0472 ± 0.0085 


0.0566 ± 0.0129 


(Tint ± Acr,-„, 


0.1325 ±0.0397 


0.1005 ±0.0261 


0.1323 ±0.0416 



metallicity dispersion seen in the four most massive globular 
clusters in M3 1 . 

In a scenario of chemical self-enrichment, the cluster metal- 
licity dispersion comes from stars originating from the gas 
trapped in the globular cluster after its formation. Therefore, one 
expects a more massive cluster to retain more gas and thus to ex- 
hibit a more prominent contribution from the second generation 
of stars. In our case, the metallicity dispersions of the three clus- 
ters are similar within the error bars. On the other hand, if their 
velocity dispersions are directly linked to their masses, G78 and 
G280 have similar masses, while G213 is less massive. With this 
information, we are not able to establish any link between the 
metallicity dispersion and the mass of the clusters. This lack of 
a correlation between the cluster metallicity dispersion and ve- 
locity dispersion could indicate that the role of the cluster mass 
in widening its RGB is not straight-forward. However it is hard 
to draw any conclusion from such a small sample of clusters. 

A second attractive explanation for the observed width of 
the RGB of these clusters would be that they are actually the 
remains of tidally stripped dwarf galaxies. Recent results on 
dwarf-globular transition objects (DGTOs -Hasegan e t al. 2005) 
and u ltra-compact dwarf (UCD) galaxies in Virgo d Jones et all 
2006) and Fornax dDrinkwater et al 2003) support the hypoth- 
esis that these objects are the remnant nuclei of threshed dE,N 
galaxies. Zinnecker et al. (1988) and Freeman (1993) have al- 
ready suggested the possibility of a> Cen being the remaining 
core of a larger entity such as a fo rmer nucleated dwarf elliptical 
galaxy. iBekki & Freemanl d2003l) suggest that the thick disk of 
our Galaxy might have been formed by the accretion of d warf 
spheroidals such as wCen. The results of M evlan et all d2001l) on 
Gl indicate this cluster could be a kind of transition step be- 
tween globular clusters and dwarf elliptical galaxies, being the 
remaining core of a dwarf galaxy whose envelope would have 
been severely pruned by tidal shocking due to the bulge and disk 
of its host g alaxy, M31. Considering that M 31 also possesses 
a thick disk dSaraiedini & Van Duvnd uOOl) and see ms to have 
suffered more tidal stri pping events than the MW (Iba ta et al 
200ll lBrown et al2 006). one would expect the presence of other 
remaining cores as "relics" of the history of interactions of this 
galaxy resembling the clusters studied in this work. In fact, both 
G213 and G280 seem to fall (at least in pro jection) on one of the 
trajectories suggested by llbata et all (|2004) for the passage of a 
satellite galaxy that might have been tidally disrupted by M3 1 . 

In order to explore this idea, we searched for any tidal fea- 
tures (tails or asymmetric envelopes) around the clusters. We 
looked for any density enhancements in the outer parts of each 
cluster to the edges of the FoV. For all three clusters there is 
no clear density excess in any particular direction that could 
indicate the presence of any tidal feature. However this visi- 
ble absence of density enhancement is not conclusive, since at 
the distance of M31, it would be difficult to detect any density 



enhancements from the star-counts of the brighter stars of the 
cluster -contrary to density enhancement detections in GCs in 
our Galaxy through star-counting methods including stars on the 
main sequence (i.e. Grillmair et al. 1995; Lee et al. 2003). 

A third possibility would be a primordial metallicity inho- 
mogeneity in the proto-cluster cloud induced by the influence of 
field stars during the formation of the globular cluster. For a long 
time the main question was whether these inhomogeneities were 
inherited at the birth of the stars that we are currently observ- 
ing (the so-called self-enrichment hypothesis that implies pol- 
lution of the intracluster gas by more massive and faster evolv- 
ing stars) or if they were generated in the course of the evo- 
lution of these objects (the so-called evolution hypothesis that 
requires non-canonical mixing inside the star itself). On the one 
hand, the classical self-enrichment scenario (Truran et al. 1991; 
Prantzos & Charbonnel 2006 -and references therein) assumes 
that stars within the GCs created the totality of the metals (light 
and heavy), and involves three stella r generation s. On the other 
hand, the pre-enrichment scenario dKraftlll994l) assumes that 
GCs were born with their current heavy metal content and in- 
vokes only two stellar generations. Both hypotheses are globally 
supported by our observations. In order to really be able to have 
conclusive arguments, we would need detailed chemical infor- 
mation (on a-, r- and s-process elements for example), while the 
color-magnitude diagnostics only give information on the cluster 
global metallicities (Z). 

In all of the above scenarios, the realm of parameter space for 
these clusters is large and we are only beginning to explore what 
the significant metallicity dispersion among their stars might be 
correlated with. In conclusion, the results presented in this work 
add to the growing body of evidence that massive globular clus- 
ters have a more complex history of star formation and chemical 
enrichment than would be expected in a traditional view of glob- 
ular clusters as simple, single-age and single-metallicity stellar 
systems. Further studies of these and similar objects will likely 
expand our understanding of galaxy assembly. 
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Fig. 8. Results of the artificial stars experiments with no statis- 
tical selection of stars for G78 (top left), G213 (top right) and 
G280 (bottom). The thin white solid line indicates the input RGB 
while the dots indicate the output (measured) RGB stars. 
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Fig. 10. Luminosity function (LF) for G78 (top), G213 (middle) 
and G280 (bottom) and their corresponding surrounding fields. 
For each cluster, the LF was derived considering the following 
inner radii, R in =5.1", 3.7", 4.2 " for G78, G213 and G280 re- 
spectively. 



Radius (arcsec) 

Fig. 9. Input and recovered magnitude difference according to 
radial position for G78 (top), G213 (middle) and G280 (bot- 
tom). For each cluster, the inner lines correspond to 5.1", 3.7" 
and 4.2" for G78, G213 and G280, respectively, while the outer 
lines correpond to 6.4", 5.3" and 5.6" for G78, G213 and 
G280,respectively. 
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Fig. 11. Luminosity function (LF) for G78 (top), G213 {middle) 
and G280 (bottom) subtracted from the contribution of their cor- 
responding surrounding fields. 
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Fig. 12. Superposi tion of the ridges and HB loci derived by 
iBrown etall d2005l) on the CMD of G78 (top), G213 (middle) 
and G280 (bottom). Thick solid line indicates the position of the 
fiducial RGB of each cluster. The straight dotted line shows the 
position of the red clump of each cluster. 
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Fig. 13. Color-magnitude diagrams (CMD) for G78 (top left), 
G213 (middle left) and G280 (bottom left) and corresponding 
surrounding fields, top right for G78, middle right for G213 and 
bottom right for G280. Grey dots indicate the stars considered 
for the fit of the RGB. Solid lines show the locus of this fiducial 
RGB. Black dots indicate the stars considered for computation of 
the RGB width. Right panels show the CMD for each surround- 
ing region. Dark dots indicate stars considered for the analysis 
of the contribution of field stars to the width of the RGB of the 
cluster. 
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Fig. 14. Color histograms for each cluster, G78 (top), G213 
(middle) and G280 (bottom), and its corresponding surrounding 
fields. Solid straight lines show the simple histogram of both 
the cluster and the surrounding field derived by simple bining 
of stars per color interval. Dots (for the cluster) and crosses (for 
the surrounding field stars) display the histograms obtained by 
weighting the contribution of each star to each color bin with its 
associated photometric error. Solid lines indicate the Gaussian 
fit for each case. 
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Fig. 15. Results of the artificial star experiments after applying 
the statistical selection of stars described in section [4] for G78 
(top left), G213 (top right) and G280 (bottom). The solid line 
indicates the input RGB while the dots indicate the output (mea- 
sured) RGB stars. Only stars falling within each cluster annulus 
were considered. 
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Fig. 16. Resulting color histogram for G78 (top), G213 (middle) 
and G280 (bottom). Filled dots display the color distribution of 
observed stars. Empty circles display the color distribution of 
recovered artificial stars. Solid and dashed lines showed the fitted 
Gaussians for each case. 
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Fig. 17. Color versus metallicit y plot for a m agnitude value of 
F814W=0.435 derived from the iBrown etal (120051) ridgelines. 
Solid line shows the fit of second order polynomial given by 
equation (1). Empty circle shows the position of the RGB of 
G78 at that magnitude value considering E(B - V)= 0.22 (top 
left panel); empty circle shows the position of the RGB of G213 
at that magnitude value considering a reddening of 0. 1 (top right 
panel); empty circle shows the position of the RGB of G280 
at that magnitude value considering a reddening of 0. 1 (bottom 
panel). For all three panels, crosses indicate the lower and upper 
limits in color and metallicity used to determine the metallicity 
dispersion of each cluster. 



